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INTRODUCTION 

The  present  work  follows  from  that  summarized  in  Ref.  1.  The 
objective  of  the  previous  work  was  to  demonstrate  by  means  of  analyses  of 
experimental  data  that  a functional  polynomial  input-output  model  for 
added  resistance  in  waves  was  a reasonable  engineering  approx imat ion. 

The  results  obtained  indicated  that  this  is  so.  It  was  demonstrated  that 
it  is  possible  to  identify  a quadratic  frequency  response  function  for 
added  resistance  from  experiments  in  both  regular  and  irregular  waves; 
and  that,  given  the  frequency  domain  response,  it  is  possible  to  synthe- 
size in  the  time  domain  at  least  the  low  frequency  nonlinear  resistance 
components . 

The  general  approach  in  the  cited  work  was  to  accept  a particular 
input-output  model  as  an  hypothesis  and  to  attempt  its  validation  by 
analyses  of  data  obtained  in  towing  tank  tests  in  irregular  waves.  The 
idea  that  a "quadratic  frequency  response  function"  for  added  resistance 
exists  and  is  meaningful  followed  from  the  hypothesized  input-output  model. 
This  idea  was  considerably  strengthened  by  utilizing  the  first  three  com- 
plete cross-bi -spectral  analyses  of  experimental  data  to  identify  the 
non-linear  responses  resulting  from  the  interactions  of  input  frequencies. 
Empirical  collapse  of  estimates  was  obtained  for  some  of  the  prominent 
features  of  the  quadratic  frequency  response  function  over  a wide  range 
of  sea  severity.  However  the  work  also  indicated  that  the  quadratic 
frequency  response  function  for  added  resistance  is  quite  complicated  — 
more  so  in  fact  than  could  be  reliably  determined  from  the  available 
experimental  data. 

Thus  in  the  cited  work  the  identification  of  the  complete  quadratic 
response  function  was  necessarily  confined  to  what  could  be  inferred  from 
existing  experimental  measurements,  and  from  the  general  (non  physical) 


1.  Dalzell,  J.F.,  "Application  of  the  Functional  Polynomial  Model  to 
the  Ship  Added  Resistance  Problem,"  Paper  presented  at  the  Eleventh 
Symposium  on  Naval  Hydrodynamics,  London,  March  1976. 
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theory.  The  weakest  part  of  the  argument  that  a meaningful  added  resist- 
ance quadratic  response  function  exists  was  that  there  existed  jto 
confirming  information  about  the  general  nature  and  behavior  of  the  complete 
function.  All  previous  theoretical  and  experimental  investigations  had 
involved  the  "mean  added  resistance  operator"  which  is  only  the  value  of 
the  quadratic  frequency  response  function  along  a very  special  line  in  the 
bi-frequency  plane.  According  to  the  experimental  results  there  appeared 
to  be  some  important  input  frequency  interaction  response  at  combinations 
of  frequencies  not  resolved  by  experiment  but  possible  in  full  scale. 

It  was  clear  that  even  an  approximate  extension  of  existing  hydro- 
dynamic  theory  for  added  resistance  to  include  the  general  response  case 
would  provide  considerable  insight  into  the  nature  of  resistance  fluctua- 
tions in  random  seas,  and,  by  providing  a physical  model,  either  lend 
credibility  to  the  previous  work  or  help  expose  its  short-comings.  This 
was  the  general  objective  of  the  present  work.  The  specific  work  carried 
out  involved  an  extension  of  existing  hydromechanic  theory,  computations 
of  quadratic  frequency  response  according  to  this  extension  for  one  of 
the  experimental  cases  of  the  previous  work,  and  a comparison  of  the 
results  with  experiment. 
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THE  INPUT-OUTPUT  MODEL 


A few  definitions  and  relations  pertaining  to  the  general  input- 
output  model  are  needed  in  the  present  work.  These  will  be  given  here 
for  convenience,  though  they  may  be  found  in  Ref.  1,  or  (in  generally 
greater  detail)  in  the  series  of  references  which  form  the  source 

material  for  Ref.  1 . 


The  general  input-output  model  is  a time  domain  model  and  is  written 
as: 

Dw(t)  = I (tl)Tl(t-t1  )dtx 

+ II  gs(t1,t8)n(t-t1)Ti(t-t2)dt1dt2  (l) 

(Omission  of  limits  on  the  integrals  signifies  limits  of  ±»  .) 

In  Eq.  (1):  0^(0  is  the  resistance  added  by  waves  (the  output) 

T) ( t ) is  a suitably  chosen  wa\  ...evation  (a  zero  mean  input) 

gx (t)  is  a linear  impulse  response 

g (tj,t  ) is  a second  degree  kernel  analogous  to  an  impulse 
response 


2.  Dalzell,  J.F.,  "Application  of  Cross-Bi-Spectral  Analysis  to  Ship 
Resistance  in  Waves,"  S IT -DL-72- 1 606,  AD749102,  Davidson  Laboratory, 
Stevens  Institute  of  Technology,  May  1972. 

3.  Dalzell,  J.F.,  "Some  Additional  Studies  of  the  Application  of  Cross- 
Bi-Spectral  Analysis  to  Ship  Resistance  in  Waves,"  S I T — D L - 7 2—  1 6^+ 1 , 
AD757363,  Davidson  Laboratory,  Stevens  Institute  of  Technology, 
December  1972. 

4.  Dalzell,  J.F.,  "Cross-B i Spectra  1 Analysis:  Application  to  Ship 

Resistance  in  Waves,"  Journal  of  Ship  Research,  Vol . 18,  No.  1, 

March  1974,  pp  62-72. 

5.  Dalzell,  J.F.,  "The  Applicability  of  the  Functional  Polynomial  Input- 
Output  Model  to  Ship  Resistance  in  Waves,"  S IT-DL-75- 1 794,  Davidson 
Laboratory,  Stevens  Institute  of  Technology,  January  1975. 
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With  respect  to  application  to  added  resistance,  the  ship  is  assumed  to 
be  proceeding  at  constant  forward  velocity  in  a wave  system  (the  input) 
which  is  defined  as  the  wave  elevation  at  a point  stationary  with  respect 
to  the  mean  position  of  the  ship;  that  is,  Eq.  (1)  is  an  encounter  domain 
model.  The  second  degree  kernel  is  assumed  to  be  symmetric  in  its  argu- 
ments, and  both  kernels  are  assumed  smooth  and  absolutely  integrable. 

The  first  term  in  Eq-.  (l)  is  the  ordinary  linear  convolution  of  a 
linear  impulse  response  function  g (tt)  with  T|(t)  • The  second  term  is  a 
double  convolution,  and  the  second  degree  kernel  ( tT , t^ ) might  be  called 
a quadratic  impulse  response.  According  to  basic  assumptions,  the  two 
impulse  responses  are  related  to  the  linear  and  quadratic  frequency 
response  functions  through  the  Fourier  transform.  These  two  sets  of 
transform  pairs  may  be  defined  as  follows: 

9,  (t)  = J e+ 1 urr  Gj  (u>)  dm 

Gi  (<")  = J e llUT  9i  (T)dT  (2) 

98<Vt8>  JJExp[+iVi  + iVs^G8  K^^d^ 

G3(ui1,aj2)  = JJ'Exp[-iulTl-iu)sT3]gs(T1,a3)dT1dT3  (3) 

The  linear  frequency  response  function,  G1(uj),  defined  by  Eq.  (2) 
is  absolutely  conventional.  The  quadratic  frequency  response  function, 

G (lD,  ,U)  ) is  defined  in  a b i -frequency  plane.  Because  the  kernel  g (t  ,t  ) 

•z  c,  2*2 

is  assumed  to  be  symmetrical  in  its  arguments: 

9s(ti>t3)  = 9s(VTi)  (4) 

and  thus: 

G8(*l»®8)  = (5) 

= GSH 

= Gs(-u)2,-uj1)  (6) 

(The  asterisk  denotes  the  complex  conjugate.) 

k 


Figure  1 summarizes  the  results  of  applying  Eq.s  (5)  and  (6)  in 
the  bi-frequency  plane  for  the  eight  possible  coordinate  positions  of  two 
frequencies  whose  absolute  values  are  a and  b.  Equation  (5)  results  in  a 
line  of  symmetry  along  the  line  0)^=1^.  Equation  (6)  results  in  a line  of 
symmetry  of  the  real  part  of  G^tu^oUg)  defined  by  uu^-uj^  (and  it  may  be 
noted  that  along  this  line  the  imaginary  part  of  the  function  is  zero). 
These  two  lines  and  the  (1^,0)  axes  divide  the  bi-frequency  plane  into 
octants,  of  which  the  two  on  either  side  of  the  positive  axis  may  be 
arbitrarily  chosen  for  reference.  The  equalities  shown  in  Figure  1 in 
the  remaining  six  octants  are  arrived  at  by  applying  Eqs.  (5)  and/or  (6). 
The  assumptions  of  symmetry  of  the  second  degree  kernel  result,  with 
Eq.  (3),  in  a complete  definition  of  G (uu,  ,uo  ) if  the  functions  are  defined 
in  any  pair  of  octants  including  a semi-axis  of  either  frequency.  Thus 
without  loss  in  generality  interpretation  of  the  quadratic  frequency 
response  needs  only  to  involve  the  octants  on  either  side  of  the  positive 
axis.  In  these  octants  uux  is  positive  and  | | a | u)g  | . 

The  interpretation  of  the  quadratic  frequency  response  function  is 
less  direct  than  for  the  linear  case,  but  can  be  approached  in  a similar 
manner.  If  in  the  linear  case  the  system  is  considered  to  be  excited  by 

71(t)  = a cosuut 

the  output  may  then  be  written: 

RejaG1  (uj)  Exp(  i u>t ) ] 

and  Gx  («)  is  interpreted  in  terms  of  normalized  amplitude  and  phase  of 
response. 

To  interpret  the  quadratic  frequency  response,  dua 1 harmonic  exci- 
tation is  necessary.  Accordingly,  it  may  be  assumed  that: 

Tl(t)  = a,  cosu)  t + a costu  t (7) 

In  accordance  with  the  previous  discussion  of  symmetry,  both  frequencies 
(u^jU)  ) are  considered  positive  and  | | s | | . The  basic  model,  Eq.  (l) 

is  good  for  any  zero-mean  excitation.  Accordingly,  Eq.  (7)  may  be  substi- 
tuted directly  in  Eq . (1).  After  some  algebra  the  final  result  for  the 
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response  to  dual  harmonic  excitation  may  be  written  as  follows: 

Ds  = Reia!  (^  ) ExpC'U)!  t)+a2G1  (uj2)  Exp(  iuost)  | 

+ |iaaG2K>-wi)+a3G3(,u^-(u2M 

+ ^ Re{a^G2(uql,oil)Exp(i2u)1t)  } 

+ \ Relaj62(u>2,m2)Exp(i2<uat)| 

+ Rejaxa2G2  (u^  ,u>2)  Exp[  i (u^+io.,)  t]} 

+ Reja1a2Ga(o)1,-aj2)Exp[i(a)1-a)2)t]{  (8) 

This  result  shows  that  the  response  of  the  quadratic  system,  Eq . (1), 
to  dual  excitation  contains,  in  general,  a shift  in  the  mean  and  components 
of  six  different  frequencies  [u^  , 2^,  2u)2,  (ut^+tUg),  and  (u^-U)  )]. 

The  first  two  terms  of  the  result  are  the  superposition  of  the  linear 
responses  at  the  excitation  frequencies.  The  third  and  fourth  terms  of 
Eq.  (8)  represent  a shift  in  the  mean.  These  terms  allow  the  identifica- 
tion of  the  mean  added  resistance  operator  as  the  value  of  G^uo^u)  ) along 
the  line  u^-u^  (or  G2  (u^  , -u^ ) ) . The  fifth  and  sixth  terms  are  the  second 
harmonic  components  (2u)a,2uJa).  Similarly,  these  terms  allow  the  identifi- 
cation of  second  harmonic  response  with  the  values  of  Ga  (u^ , ) along  the 

1 i ne  U3s=u\  (or  G2(u)1,uj1))  . 

The  seventh  and  eighth  terms  of  Eq . (8)  pertain  to  the  bi-frequency 
plane  in  general.  The  seventh  term  is  the  response  at  frequency  (u^+u^)  ; 
that  is,  Gg  (0^,111  ) expresses  the  normalized  response  in  the  sum  frequency 
due  to  non-linear  interactions.  Similarly,  the  eighth  term  involves 
response  at  frequency  (u^-u^);  that  is,  G^u^,-^)  is  the  normalized 
response  in  the  difference  frequency.  In  terms  of  Figure  1,  the  octant  of 
the  b i -frequency  plane  above  the  positive  u>  axis  corresponds  to  the  por- 
tion of  the  quadratic  frequency  response  function  which  defines  sum- 
frequency  interactions,  and  the  octant  below  the  positive  axis  corres- 
ponds to^the  portion  of  the  function  which  defines  difference  frequency 
i nteract  Tons . 
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Equations  (7)  and  (8)  together  afforded  a direct  way  to  proceed 
in  the  present  investigation.  In  particular,  what  is  required  for  a 
hydromechanical  solution  for  the  quadratic  frequency  response  is  to 
solve  the  hydromechanical  problem  for  dual  regular  wave  excitation. 
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GENERAL  HYDRODYNAMIC  FORMULATION 


The  initial  part  of  the  theoretical  formulation  is  due  to  Salvesen 
In  that  reference  the  hydrodynamic  force  on  a body  in  the  free  surface  is 


given  as: 


f = -p  Ht  u * Kds  -p  n RS v$-  i i v$  ?*}ds 


SB+SF 


In  the  above: 


F = hydrodynamic  force 
$ = total  velocity  potential 

n = outward  unit  normal  vector 
p = mass  dens i ty 
t = t ime 

Sg  = submerged  surface  of  the  body 

Sp  = portion  of  the  free  surface  inside  a far  field 

control  surface,  S 
7 00 

No  approximations  were  made  in  deriving  Eq.  (9).  The  expression  is  taken 
as  being  exact  within  potential  flow  theory. 

If  it  is  assumed  that  the  potential  $ is  composed  of  a steady  state 
wave  resistance  potential,  ?^,and  a time  dependent  potential  due  to  waves, 

*T  ' 

$ = $s  + $T  (10) 


6.  Salvesen,  N.,  "Second-Order  Steady-State  Forces  and  Moments  on  Surface 
Ships  in  Oblique  Regular  Waves,"  Paper  22,  International  Symposium  on 
the  Dynamics  of  Marine  Vehicles  and  Structures  in  Waves,  London,  April 
1974. 
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Then  Eq.  (9)  becomes: 


-P 


-P 


-P 


7t  If  $T  " 

r _i_r  1 


dS 


SB+SF 


{•f  [jT  °‘t  + ST  ’‘s'  (,*s ' «T)s]ds 


.9$. 


If  r^T  V$T  - ~ | VfT!2n]dS 


00 


In  the  present  case  the  interest  is  in  forces  which  are  related  to 
waves.  Accordingly,  the  first  integral  in  Eq.  (11)  may  be  ignored  as 
being  constant  and  independent  of  waves. 

If  it  is  assumed  that  the  time  dependent  potential  is  linear  with 
respect  to  waves,  and  is  expressible  at  most  in  the  form  of  sums  of  con- 
tributions from  wave  components  of  different  frequencies,  the  second  and 
third  integrals  make  contributions  to  force  which  are  also  linear  with 
waves.  Under  this  assumption  these  integrals  are  also  not  of  interest  in 
the  present  problem. 

The  last  integral  in  Eq.  (11)  involves  products  of  time  dependent 
potentials.  It  is  thus  inherently  a quadratic  nonlinearity  even  under 
the  assumption  of  linear  potentials.  The  cardinal  assumption  in  the 
present  work,  and  that  of  Salvesen  , is  that  the  dominant  contribution 
to  the  quadratic  or  second  order  non-linearities  comes  from  the  combina- 
tion of  linear  or  first  order  potentials  in  the  fourth  integral  of  Eq.  (11) 
and  not  from  nonlinearities  of  the  potentials  themselves. 

The  quadratic  part  of  the  added  resistance  or  drifting  force  is 
thus  taken  to  be: 


where  is  the  first  order  time  dependent  potential.  This  expression  is 
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essentially  the  same  as  that  taken  by  Salvensen  in  formulating  the  mean 
drifting  force,  and  his  development  can  be  followed  further. 

Writing  the  time  dependent  potential  as  the  sum  of  an  incident 
wave  potential,  and  a potential,  due  to  the  body  disturbance, 
including  diffraction  effects: 


*T  = $B  + *|  03) 

A result  attributed  to  J.N.  Newman  is  used  by  Salvesen  to  put  Eq.  (12) 
in  the  following  form: 


R = pJT 

Sco 


* 9 

L$B  Bn  " 


fV] 

3n  _ 


(V*,  + £ v§B)dS 


(14) 


Salvesen  made  the  assumption  that  the  body  is  a weak  scatterer 
so  that:  This  assumption  is  considered  to  be  questionable 

for  oblique  ship-wave  headings.  However,  the  assumption  is  considered 
quite  reasonable  for  head  seas  (the  case  of  present  interest)  and 
Salvesen's  calculations  bear  this  out.  In  order  to  reduce  the  complexity 
of  the  present  analysis,  the  weak  scattering  assumption  was  accepted. 
Accordingly  Eq.  (14)  becomes: 


* ■ K b ■ *r]  ,4I  dS 


This  equation  corresponds  to  Eq.  (32)  of  Salves 


en 


05) 


Application  of  Green's  Theorem  changes  the  integration  in  Eq.  (15) 
to  one  over  the  body  surface: 


« ■ -oil  [*.  i 


^B~l 
Bn  J I 


dS 


(16) 


Now  the  present  approach  to  the  interpretation  of  the  quadratic 
response  function  involves  the  computation  of  response  to  two  harmonic 
waves  of  different  frequencies.  Given  the  assumption  of  linear  poten- 
tials, this  involves  a simple  superposition  of  potentials,  which  does  not 
upset  the  conclusions  made  in  the  discussion  of  Eq.  (11).  Accordingly, 
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the  time  dependent  incident  and  body  potentials  corresponding  to  two 
incident  waves  of  frequencies  tu.  and  ojj  were  assumed  in  the  following 
form: 

$5  = Re  ^cpB(w.)  Exp[-iu).t]  + Exp[- iuh  t]~j 

= \ [cPB(u,i)ExP["'“}it]  + ^(uij)  Exp[  iu>.t] 

j.  — 

+ cpB((Uj.)  Exp[-iujjt]  + 9g(u>j)  Expfiuijt]  (17) 

$1  = Re  ^cp,  (uj.)  Exp[-iu)jt]  + cp|(u,j)Exp[-iU)jt]] 

= \ 1*9,  (“>;)  ExpC-ioi.t]  + cp“(a».)  Exp[  iiu.t] 

+ cpj  ((Uj)  Exp[-iio.t]  + ^(uij)  Exp[  iu)jt]j  (18) 

The  complex  potential  functions  cp ^ (u^)  Exp|"- i u^t]  and  cp^(u^)  Exp[-iu)^t]  have 

the  conventional  interpretation  in  ship  motions  theory,  they  correspond  to 

t h 

the  incident  and  body  potentials  in  the  k regular  wave  train,  the  ampli- 
tude of  the  wave  train  is  imbedded  in  the  cp  and  cp  functions,  and  the 

I D 

frequency  is  encounter  frequency — not  wave  frequency. 

Substituting  Eqs . (17)  and  (18)  into  Eq.  (l6)  and  gathering  terms 
there  results: 

\ = Re  [Hj  (u).,-oo.)  + Hx  (uJj  , -UJj ) 

+ Re  !"ht  (u)j,u>.)  Exp[-i2u).t]J 
+ Re  (to. , iu j ) Exp[- i2tu.  t]  j 

+ Re  ( UJ j , UJj ) + (ojj ,(«.)  j-Exp[- i t (m.+U)j ) ] j 

+ Re  (u).,-u>j)  + H"(u)j,-u).)jExp[-it(u).-u)J.)]j  (19) 

where: 

R^  = force  due  to  two  superimposed  regular  wave  trains 
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= • 2 M fp{CpB(u,i)  h " h cPB(u,i)^vCPl(u,j)]dS  (20) 

SB  U 


V- 


c 


Now  forming  a new  function,  H (w.,U).): 

s 1 J 


H3(«>i,UJj)  « Hx  (<*>.,<«.)  + (lUj  , u). ) 

Hs  (<";.-«•)  = H4  (»,,-«.)  + H*  (ujj  , -U).  ) 


(22) 

(23) 


This  function  has  the  same  symmetries  with  respect  to  frequency  arguments 
as  the  quadratic  frequency  response  function,  ).  Specifically, 

as  i n Eqs . (5)  and  (6) : 

H,  (».,»,)  = H2(u).,U)j) 
and,  since  cpB(-iu)  = cp“ (uj)  : 

Ha  (*®i  ,Uj  ) = Hs  (-u,i  ^ _0Jj  ) 

= H8(-(Uj.,-U).) 

Substituting  Eqs.  (22)  and  (23)  in  Eq.  (19): 

I 


Rs  = 2 [Hs(u,i'’tui)  + MVV] 


+ \ Re  [” Hs  ( u> . ,(U.)  Exp[-i2ui.t]l 


+ \ Re  | Hs  (uj.,u).)  Exp[-i2ui.t]j 
+ Re  [" H2  (tu. , m. ) Exp[-  i t (u).+u!j ) 

+ Re  [H8(u).,-u)j)Exp[-it(u).-U)j)]j 


(24) 


Thus  the  second  order  force  induced  by  two  regular  wave  trains  has 
been  brought  into  essentially  the  same  form  as  the  last  six  terms  of  Eq.  (8) 
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These  last  six  terms  involve  the  contribution  of  the  quadratic  non-linearity 
in  the  response  of  the  general  input-output  model  to  dual  harmonic  excita- 
tion. 

Noting  from  Eq.  (23)  that  Hs(id.,-iu.)  is  purely  real,  a comparison 
of  Eqs.  (24)  and  (8)  yields: 

Ga(u>,,-ui,)»  —i"  Hs(uJ.,-u>i)  (25) 


where  a.  is  the  amplitude  of  wave  component  of  frequency  a).  . 
Continuing  in  the  same  way: 


Gs(oui^i)  = IJ  H3,(("i^i)  (26) 

= TT:  Ha(«*»i>»|)  (27)  . 

and:  , * 

Ga(uji>-u,j)  Hs(tui»-“)j)  (28) 

J 1 j J 

Thus  the  formulation  of  the  computation  of  the  quadratic  frequency 
response  function  is  complete.  For  each  pair  of  regular  waves  of  ampli- 
tudes a.  and  a^ , and  frequencies  u).  and  to.,  Eqs.  (20)  and  (21)  must  be 
evaluated  to  form  ^ (u).  ,0)^. ) , Ht  (uh,U).),  H1(u).,-u)j.)  and  H1  (uj.  , -u>. ) . From 
these  results  Eqs.  (22),  (23),  (27)  and  (28)  may  be  evaluated  in  turn  to 
yield  the  desired  estimates. 

One  observation  is  in  order  before  continuing.  The  present  develop- 
ment involves  relatively  little  specialization  to  added  resistance. 
Computation  of  quadratic  frequency  response  functions  may  be  made  accord- 
ing to  the  above  recipe  for  other  forces  within  the  assumption  of  linearized 
potential  flow  for  a body  which  may  be  considered  a weak  scatterer.  Removal 
of  the  weak  scattering  assumption  appears  feasible  at  the  expense  of  a 
moderate  increase  in  complexity  of  the  resulting  hydrodynamic  expressions. 

Finally,  it  is  of  interest  to  compare  the  present  result  with  the 
corresponding  stage  in  Salvesen's  development  for  the  mean  force  induced 
by  a single  wave.  To  do  this  let  the  amplitude  of  the  j**1  wave  component 
go  to  zero.  Since  the  amplitude  is  imbedded  in  the  potentials: 


O 


I 
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= ^((uJj)  = 0 and  thus  only  the  first  and  third  terms  remain  in 
Eq.  (19).  Taking  the  time  mean  of  the  first  and  third  terms  of  Eq.  (19) 
leaves  only  the  first  term: 

(u>.,-u).)  ] 

Apart  from  the  explicit  notation  that  the  real  part  is  to  be  taken,  and 
some  purely  notational  differences,  this  result  is  exactly  the  same  as 
Salvesen's  Eq.  (3*t)  for  the  mean  force. 
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ANALYTICAL  EVALUATION  OF  THE  QUADRATIC 
FREQUENCY  RESPONSE  FUNCTION  FOR  ADDED  RESISTANCE 


While  there  is  no  fundamental  objection  to  the  development  of  a 
vector  quadratic  frequency  response  function  corresponding  to  the  vector 
force  represented  by  the  equations  of  the  last  section,  the  objectives 
of  the  present  work  involved  only  a comparison  with  experimental  results 
for  added  resistance  in  head  seas,  and  it  was  accordingly  a worth-while 
simplification  to  consider  only  the  longitudinal  component  of  force  for 
the  head  sea  case. 

In  general,  the  evaluation  of  the  functions,  Eqs.  (20)  and  (21), 
was  carried  out  in  a manner  similar  to  that  employed  by  Salvesen°  in  that 
a strip  method  was  used,  and  the  same  or  very  similar  hydrodynamic  assump- 
tions were  employed.  Some  differences  do  exist  because  the  strip  method 
used  was  that  of  Ref.  7. 

The  most  important  difference  between  this  method  and  others 

(Ref.  6 for  example)  is  in  the  evaluation  of  the  diffraction  induced  wave 

exciting  forces  and  moments.  The  method  of  Ref.  7 is  essentially  an 

extension  to  forward  speed  of  the  method  for  zero  speed  developed  in  Ref.  8. 

In  both,  the  diffraction  induced  forces  and  moments  as  well  as  the  motion 

g 

induced  forces  and  moments  are  evaluated  by  the  Frank  close  fit  method. 
Because  the  work  of  Ref.  7 has  not  as  yet  received  wide  circulation,  a 
summary  of  the  method  is  included  herein  as  Appendix  A. 


7.  Kim,  C.H.,  "Calculation  of  Motions  and  Loads  of  a Ship  Uniformly 
Advancing  in  Oblique  Regular  Waves,"  Technical  Memorandum  SIT-DL-7^- 
1 66,  March  197*+. 

8.  Kim,  C.H.  and  Chou,  F.,  "Wave-Exciting  Forces  and  Moments  on  an  Ocean 
Platform  in  Oblique  Seas,"  Paper  No.  OTC  1180,  Offshore  Technology 
Conference,  April  1970.  Also  SIT-0E-71-Report  6,  September  1971. 

9.  Frank,  W.,  "On  the  Oscillation  of  Cylinders  in  or  Below  the  Free 
Surface  of  Deep  Fluids,"  NSRDC  Report  2375>  October  1967. 
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Figure  2 indicates  the  usual  right-handed  coordinate  system  used 
in  the  present  development.  It  is  described  in  more  detail  in  Appendix  A. 
In  the  present  case  the  ship  is  presumed  to  be  restrained  in  surge  and 
to  be  proceeding  at  constant  velocity  U in  two  superimposed  head  waves 
(p,  = rr  in  Figure  2) . 


In  a manner  consistent  with  the  form  of  the  incident  wave  potential 
assumed  in  Eq.  (18),  the  two  regular  waves  are  represented  by  the  real 
parts  of: 

-i(v  x+au.t) 

o.  i 


and : 


C,  = 


a.  e 
J 


-i  (v  x+u).  t) 

°j  J 


(29) 


Here  the  i,j  subscripts  indicate  two  different  waves,  and 

a.  = wave  amplitude 

v = a)  ./g,  the  wave  number 
o . o i ' 

tUj  = u)Q  + Uvq  , frequency  of  encounter 


The  incident  wave  potential,  9((u)j),  for  the  jth  head  wave  trjin  is 
given  by  (omitting  the  time  factor  in  accordance  with  Eq.  18): 


ga.  V(z“ix) 
9.  (oj.)  = -i  — L e J 

1 J “o. 

J 


(30: 


i.  L 

The  gradient  of  the  incident  potential  for  the  j wave  is 


V9,  (u>j)  = [T(-ivo  ) + kvQ  ]cp.(u).)  (3 1 ; 

j j J 

In  order  to  specialize  the  analysis  to  added  resistance,  only  the 
longitudinal  component  of  Eq.  (31)  is  taken.  Accordingly,  in  the  evalua- 
tion of  Eqs . (20)  and  (21)  the  gradients  of  the  incident  wave  potentials 
may  be  replaced  as  follows: 


yep,  (u>  ) - -ivQ  cp,  (uj  ) 

j J 


v9*(«a>j)  - ivn  cp* (cu. ) 


(32' 
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with  the  understanding  that  the  result  will  apply  only  to  the  longitudi- 
nal component  of  force. 

The  body  potential  for  the  i^  wave  train  is  taken  as  the  sum  of 

the  radiation  and  diffraction  potentials.  Because  of  the  heading,  only 

vertical  modes  of  motion  are  considered.  Thus  for  the  response  to  the 
. th 

i wave  tram: 


V<u  i)  = 2 V(m)(u>.)9(m)(u).)  + a.cpn(,o.)e 

m=3,5  ' ' U ' 


■' V* 


where: 


V m (uk)  = velocity  amplitude  of  heave  (m=3)  and  pitch  (m=5) 

cp^(u).)  = radiation  potential  due  to  heaving  and  pitching 
motions,  each  per  unit  velocity  amplitude 

<pp(<«j)  = the  diffraction  potential  per  unit  wave  amplitude 
in  the  i wave  train  (see  Appendix  A) 

V (3)  (u)j ) = - iu).£  (u). ) 

(U3.)  = - iu>.  \|f  (uj.)  (3J 

and  C and  i|r  indicate  the  heave  and  pitch  displacement,  respectively. 

t h 

The  diffraction  potential  for  the  i wave  train,  cpQ,  (Appendix  A) 
is  evaluated  by  satisfying  the  kinematical  boundary  condition: 


[aiCPD(u)i) 


- IV  x 

°i  1 


_ = - ^ 9,(0).)  on  SE 


Substititing  Eqs.  (32),  (33),  and  (35)  into  Eqs . (20)  and  (21) 


there  results: 


6 


i 
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H,  = 1 Vo p{l  dS 

J 5d 


-IV„  x 


+ ai  n cpD(u,i)< 


o,x  39,  (U».) 


, \ a<p^ (w.) 

Vf,5  V ” (“i>  l-TsS i-  dP.fuKjdS 


r 


q>!  (u)j)dS^ 


«,(«■-.)  =i,  ,{-  z »«•>(.,)  JT  <p(”)(«,) 

J j 1 m=3,5  s 


Scpj'(UJ  ) 


„ "VX  3<P*(UJ.) 

- ai  II  1 ^ — L ds 


’!  jH  i ‘ 

B 


+ S V(m)(*.)  JJ=- 
m=3,5  1 * Sn 


acp^m^  (ou.) 


- 9 ,”  ( tu  j ) d S 


Sep  (uj  ) 

<-J>d*}  I 

B 

In  the  radiation  problem  the  kinematical  boundary  conditions  for 
heaving  and  pitching  are  given  by 

s?(3)  _ 

■r - n 

an  - 

f? -\n(x  + J^)  on  SB  ( 


where  ri  is  the  z-component  of  the  unit  normal  on  the  hull  surface  (in 
fact,  it  is  on  a section  surface  in  the  stripwise  calculation)  pointing 
into  the  water. 


r 
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From  the  above  we  have  the  relationship  between  cpu'  and  cp^  , 


(u>  ) = -(x  + (u>.) 


The  partial  of  the  incident  wave  potential  is  evaluated  as: 


dcp,  (u>.) 


- mo  cp,  Co,  ) 

J ‘ 


where  in  the  slender  body  approximation  the  x-component  of  the  unit  normal 


is  zero. 


Further,  from  Eq.  (38): 


89^  («J.)  , s.  , . 

«,  T7 - «•>>  * iUm  m 

i an  i 


n<3>-n,  b(3).0<  m(5),„ 


accord ing  to  the  variant  of  Stokes'  theorem,  we  have: 


where 


JJV 


8cp(m)  (u>.) 


9,(w.)ds  = JJ  fou  n(m)  + iUm(m))9,  (ui.)dS  (43) 
1 J C 1 I J 


It  may  be  seen  that  substitution  of  Eq.  (40)  through  (43)  in 
Eqs.  (36)  and  (37)  results  in  a factor  of  n in  each  term.  Accordingly 
the  double  integral  over  Sg  in  the  equations  may  be  replaced  by  a strip- 
wise  integration: 

JJ  H8S  = - J dx  f dy  (L 

SB  L C(x) 

where  C(x)  = a transverse  section  contour 
y = a transverse  coordinate 
L = hul 1 length 

Further  transformations  of  the  terms  in  Eqs.  (36)  and  (37)  are 
effected  as  follows:  The  sectional  heave  added  mass  m"(u).)  and  damping 

coefficient  N(uj.)  are  obtained  from  the  following  relation. 


C 
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pw.  J 9 (uj.)dy  = u).m"(u).)  + i N (co. ) (45 

C(x)  ...  . 

The  sectional  heave  diffraction  force  per  unit  wave  amplitude  (see 
Appendix  A)  is  given  by: 


ipuj.  r 

' Ctx) 


cP0(u,i)dy  = fHu(«u.) 


The  sectional  restoring  force  per  unit  heave  displacement  is  given 


pg  j dy  = pgB  (A 

C (x) 

where  B = beam  of  the  section  in  the  load  waterplane. 

By  taking  the  average  section  draft  T at  x in  the  integral  of  the 
exponential,  an  approximate  expression  can  be  derived: 


v z -v  T 

o.  o. 

f e J dy  = e J , dy  (4£ 

C(x)  C (x) 

, = wetted  sectional  area 

where  T = — — : r 

wetted  section  beam 

Upon  substitution  of  Eqs.  (30),  (34),  (38)  through  (40),  and  (43)  through 
(48)  into  Eqs.  (36)  and  (37)  there  resulted  expressions  which  could  be 
evaluated  with  the  aid  of  existing  ship  motion  programs. ^ 
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fCK) 


-iv  x -v  T 
o.  o. 


1 r^i'  . u;  u; 

2 VjW  FIT-  J e J e J [u>.m"(u).)  + «N(ai)]d> 


J J 


- i v x -v  T 


'K'Uj)  y O.  'O. 

~ — J (x+  TaT^6  J e J Cu>.m"((u.)+  iN(iu.)]dx>- 


i L 


‘V  “i(vo.+vo>  -Vof 

* I Vjvo  F if  * ' J e J fHD(«.)dx 


o a.a.v  t-1  + — L 

2 .jo.  L v 


J 


f C(*».)  -ivo.X  -o.T 

~ — J e J e J pgB  dx 


i L 


*(«*.)  -ivo.*  ‘VT 


a . 


j xe  J e J pgB  dx} 


i L 


“o.“o.  •i(V*''o>  -<''o+''0)T 

♦ivj-t-1;*  ' J • ’ j 


pgB  dx  (49) 
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Hi 


iv  x -v  T 
o.  o. 


t;  aia:v0  “o  1 a"‘"  I e j e j [uj.m" (UJ. ) + iN(ou.)]dx 

J j j i L 


[ / + JL  'V°iX  e °j  tu,im"(u)i)+  'N(u).)]dxj 
J v iu). ) e J 


-v  T 


a . 

1 


<v 


o. 


- i (v  -V  )x  -v  T 

o.  0/  o.  n 

+ — a.a.v  — L f e ' J e ^ fu  (<D.)dx 

2 1 j o . u).  J H 1 

J • L 


oj.-uj.-, 

_J J. 


+ 7r  a.a.v  I 1 + 

2 1 j o.L  U) 

J O; 


f C(u,i}  r 

r 17“ 


iv  x -v  T 

o.  O. 

! 1 e pgB  dx 


i v x -v  T 
o.  o. 


i(u>.) 

Y I p I I 

+ J xe  e J pgB  dx* 


+ 2 aiaj 


U)  u> 
o.  o. 

! L 


•i  (v 


I e 

L 


-v  s 
o.  o.)x 


-(v  4v  )f 

o.  o. 

e ' J pgB  dx  (50) 


For  present  purposes  the  expressions  (Eqs.  49  and  50)  for  ^(ui.jio.) 
and  ^ (uj.  , ~(Uj ) may  be  used  in  place  of  Eqs.  (20)  and  (21)  in  the  expres- 
sions for  Hg  (uu. , u). ) , etc.,  Eqs.  (22)  and  (23),  with  the  understanding  that 
the  resul t relates  only  to  the  quadratic  frequency  response  function  for 
the  longitudinal  component  of  force.  Equations  (49)  and  (50)  contain  the 


product  of  the  regular  wave  amplitudes  (a;aj)  as  a factor.  Thus  the 
functions  Ha(u).,u)j),  Eqs.  (22)  and  (23),  have  the  same  factor,  and  when 
the  quadratic  frequency  response  function  is  evaluated  by  Eqs.  (25) 


through  (28)  the  factor  cancels  — meaning  that  Eqs.  (49)  and  (50)  may 
be  evaluated  for  unit  wave  amplitudes. 


The  dimension  of  the  quadratic  frequency  response  function  is 
( force/ 1 eng 1 7) . In  the  reduction  of  experimental  data  (Ref.  1-5)  length 

was  non-d i mens ional i zed  with  model  length,  L,  forces  were  non-d imens i on- 
alized  with  model  displacement.  A,  resistance  was  considered  to  be  a 
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positive  force,  and  frequency  was  non-dimens ional ized  wi th  /2rrg/L  . 

Thus  the  last  details  in  generating  theoretical  estimates  comparable  to 
those  in  Refs.  1-5  are  to  change  the  sign  of  the  results  of  Eqs . (25) 
through  (28)  and  multiply  by  L2/A,  and  compensate  for  the  frequency  non- 
d imens iona 1 iz i ng  constant. 

As  may  be  noted  from  Eqs.  (22)  and  (23),  the  answer  sought  for 
bi-frequency  u).,<u.  involves- the  evaluation  of  Eqs.  (49)  and  (50)  as 
written,  and  with  the  i j subscripts  interchanged.  In  addition,  it  was 
required  to  develop  the  theoretical  estimates  over  a region  in  the  tu  ,to 
plane,  Figure  1.  This  means  that  results  from  the  basic  ship  motions 
computation  for  each  frequency  must  be  used  over  and  over  again.  Accord- 
ingly the  actual  computation  was  divided  into  two  steps.  Before  commenc- 
ing the  first  step  the  required  extent  and  fineness  of  coverage  of  the 
bi-frequency  plane  are  decided  upon,  and  a value  of  6cu  and  K are  chosen. 
These  parameters  define  the  discrete  points  in  the  plane  at  which  esti- 
mates of  the  quadratic  frequency  response  will  be  generated;  that  is, 
estimates  are  thus  specified  to  be  made  for  all  combinations  of: 


<^=0, 6u),  • • • 'K6uj 


u>2  = -K6u),  • • *,-6uj, 0,6ur  • "KSo) 


Given  a forward  ship  speed,  U,  the  wave  numbers  corresponding  to  each  of 
the  (K+l)  encounter  frequencies  may  be  solved  for. 

In  the  first  step  of  the  computation  the  ship  motions  problem  is 
solved  for  each  of  the  (K+l)  discrete  encounter  frequencies.  Essentially, 
the  sectional  added  mass,  damping,  restoring  force  and  exciting  force  are 
developed  for  each  strip  and  these  are  integrated  to  result  in  estimates 
of  pitch  and  heave  amplitudes.  All  of  these  results  are  stored  so  that 
the  result  of  the  first  computat iona 1 step  is  a "library"  of  sectional 
added  mass,  damping,  etc.,  for  all  the  specified  encounter  frequencies. 
Results  for  "zero"  frequency  are  approximated  by  those  for  a very  low 
frequency. 

In  the  second  step  of  the  computation  Eqs.  (49)  and  (50)  are 
systematically  evaluated  for  each  combination  of  the  discrete  frequencies 
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by  "looking  up"  the  previously  computed  data  from  the  first  computational 
step,  and  the  estimation  of  the  quadratic  frequency  response  continues  in 

accordance  with  Eqs.  (22),  (23)  and  (25)  through  (28). 
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COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  ESTIMATES 
OF  THE  MEAN  ADDED  RESISTANCE  OPERATOR 

As  previously  noted,  the  mean  added  resistance  operator  is  the 
value  of  the  quadratic  frequency  response  function  along  the  line  (^=-<1^, 
and  it  is  this  special  section  of  the  function  for  which  there  is  a great 
deal  of  confirming  information.  Thus  the  initial  comparisons  with  experi- 
ment of  the  present  analytical  method  were  of  the  mean  operator.  The  model 

used  in  the  experiments  of  Ref.  1-5  was  of  the  Series  60  0.60  C parent, 

B 

and  this  model  was  also  one  of  those  used  in  experiments  reported  by 
1 0* 

Strom-Tejsen,  et  al 

Figure  3 shows  results  from  the  present  analysis  plotted  on  what 
is  effectively  Figure  8a  of  Ref.  10.  The  work  in  Ref.  10  on  the  Series  60, 
0.60  block  model  had  involved  two  speeds  (F roude  numbers  0.266  and  0.283). 
The  evaluation  of  the  present  analytical  method  was  carried  out  at  the 
mean  of  the  two  speeds,  Fn  = 0.274,  and  the  results  were  non-d i mens ional i zed 
in  accordance  with  the  conventions  employed  in  Ref.  10.  The  present  results 
are  clearly  closer  to  experiment  for  this  case  than  are  the  results  of  the 
analytical  methods  of  Ref.  10.  In  addition,  at  high  encounter  frequencies 
the  present  method  indicates  the  leveling  out  trend  often  observed  in 
experimental  added  resistance  operators. 

The  experimental  work  in  Refs.  1-5  had  also  involved  two  speeds, 

F =0.15  and  0.20.  The  lower  of  the  two  speeds  was  selected  and  the 
quadratic  frequency  response  function  was  evaluated.  For  this  case  Fig- 
ure 4 indicates  the  comparison  with  experiment  of  the  mean  added  resistance 
operator  from  the  present  analytical  method.  Except  for  the  addition  of 
the  analytical  result,  Figure  4 is  also  Figure  4 of  Ref.  1 or  Figure  2 of 
Ref.  5.  In  this  figure  the  oe  is  a non-dimensional  encounter  frequency: 

oe  = U)e//2rTg7L 

where  ui  is  encounter  frequency  in  radians/second. 

6 

10.  Strom-Tejsen,  J.,  Yeh,  H.Y.H.,  and  Moran,  D.D.,  "Added  Resistance  in 
Waves,"  SNAME  Vol.  81,  pp  109,  1973. 
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The  added  resistance  operator  is  non-dimensional i zed  by  L /A  as  noted  in 
the  previous  section.  Two  kinds  of  experimental  data  are  included  in  the 
figure.  The  regular  wave  envelope  indicates  the  magnitude  and  scatter  of 
the  results  of  regular  wave  experiments  reported  in  Ref.  2.  The  lines 
broken  by  various  symbols  are  the  results  of  cross-bi-spectral  analyses 
of  data  obtained  in  three  irregular  seas.  If  the  ship  is  assumed  to  be 
500  feet  in  length,  "Sea  State  A"  corresponds  to  a significant  wave  height 
of  9 feet,  "Sea  State  B"  corresponds  to  a significant  height  of  18  feet, 
and  "Sea  State  C"  corresponds  to  a significant  height  of  33  feet. 

In  Figure  4 the  agreement  with  experiment  of  the  analytical  esti- 
mates of  the  mean  added  resistance  operator  appears  to  be  somewhat  better 
than  that  shown  in  Figure  3.  In  this  case  most  of  the  analytical  result 
is  within  the  regular  wave  envelope. 

Figures  3 and  4 clearly  demonstrate  the  validity  of  the  present 
analytical  method  relative  to  experimental  data  obtained  both  at  Davidson 
Laboratory  and  elsewhere.  While  it  was  not  the  purpose  of  the  present 
work  to  refine  estimating  methods  for  the  mean  added  resistance  operator, 
the  present  method,  which  involves  a different  approach  to  computation  of 
exciting  forces,  appears  to  have  promise  in  this  direction.  For  the 
particular  ship  form  and  speeds  shown,  the  differences  between  the  present 
method  and  the  mean  experimental  result  appear  to  be  of  the  same  order  as 
the  scatter  in  the  experiments  (which  are  noted  for  their  difficulty). 
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COMPARISON  OF  ANALYTICAL  AND  EXPERIMENTAL  ESTIMATES 
OF  THE  QUADRATIC  FREQUENCY  RESPONSE  FUNCTION 


In  order  to  compare  the  analytical  estimates  of  the  complete  quad- 
ratic frequency  response  function  with  those  from  experiment  it  is 
advantageous  to  introduce  the  frequency  transformation  of  Ref.  1.  This 
transformation  maps  the  plane  (Fig.l)  into  a sum  and  difference 

frequency  plane  in  accordance  with: 

ftj  = = difference  frequency 

= uij+u>s  = sum  frequency  (51 

TheQ1^Q3  plane  is  shown  mapped  into  the  0^,0^  plane  in  Figure  5.  The 
and  Qa  axes  are  coincident  with  the  lines  of  symmetry  noted  in  the 
discussion  of  Figure  1.  In  accordance  with  the  discussion  of  Figure  1 
the  fundamental  symmetry  of  the  quadratic  frequency  response  function 
allows  complete  definition  if  the  function  is  known  in  a quadrant  of  the 
plane  defined  by  the  lines  of  symmetry.  The  quadrant  of  the  ,Q  plane 
(fig. 5)  in  which  both  sum  and  difference  frequencies  are  positive  corres- 
ponds to  the  quadrant  selected  earlier  in  interpreting  the  response  to 
dual  harmonic  excitation.  This  quadrant  will  be  used  in  the  comparsions. 
Also  in  accordance  with  the  earlier  discussion  of  dual  harmonic  excita- 
tion, the  sum  frequency,  corresponds, and  is  numerically  equal  to, the 

output  frequency.  Any  combination  of  input  frequencies  id  and  uj  which 

1 s 

are  on  the  lineQ^  = constant  will  generate  a component  at  that  frequency. 
For  example,  the  mean  added  resistance  operator  involves  output  frequency 
of  zero  and  is  located  on  the  axis. 

For  comparison  purposes  the  frequency  transformation  involves  no 
changes  to  the  magnitude  of  the  analytical  estimates  of  the  quadratic 
frequency  response  function,  merely  a somewhat  more  convenient  way  of 
looking  at  the  function.  The  analytical  results  for  the  Series  60 
0.60  block  model  at  Froude  Number  of  0.15  were  non-d i mens iona 1 i zed  in 
accordance  with  the  experimental  convention  noted  in  a previous  section. 
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The  discrete  set  of  frequencies  at  which  the  function  was  evaluated  were 
chosen  so  that  results  would  be  available  for  the  exact  non-dimensional 
frequencies  used  in  the  experimental  identification  process.  The  non- 
dimensional  frequencies  are  denoted  in  Ref.  1 by  a and  instead  of 
and  u)2  . Since  the  ct's  differ  from  the  iu's  only  by  a scale  factor 
ol  and  as  may  be  substituted  for  U)T  and  in  Eq.  (50  and  in  fact  in  all 

previous  equations.  For  present  purposes,  however,  and  ^ will  be  taken 

to  denote  non-dimensional  frequency. 


Before  proceeding  with  detailed  comparisons  an  overview  of  the 
magnitude  of  the  analytical  result  is  instructive.  Figure  6 is  an  iso- 
metric plot  of  the  modulus  of  the  analytically  determined  quadratic 
frequency  response  function.  The  "view"  is  toward  the  origin  of  the  bi- 
frequency plane  from  a position  above  the  axis,  Figure  5.  (The  differ- 
ence frequency  axis  (0  ) is  to  the  left,  the  sum  frequency  axis  is  to  the 
right,  and  the  modulus  of  the  function  is  "vertical".)  The  lines  defining 
the  function  contours  are  in  the  nature  of  section  lines  which  are  parallel 
to  the  Qx  and  axes.  Each  intersection  corresponds  to  an  analytically 
computed  point,  and  the  connecting  line  segments  are  straight.  There  were 
actually  about  twice  as  many  points  computed  as  are  shown.  It  is  a con- 
sequence of  the  point  spacing  chosen  for  plotting  that  the  function  appears 
to  be  irregular  near  maxima.  Inspection  of  the  intermediate  numerical 
results  indicated  that  the  computation  defines  a smooth  surface.  In  the 
figure  the  plotted  domain  of  the  bi-frequency  plane  appears  triangular. 

This  is  the  result  of  computing  the  function  to  a maximum  absolute  input 
frequency  of  2.62.  The  r ight-to- 1 ef t section  at  the  front  of  the  plot  is 
shown  in  bolder  lines  with  vertical  lines  defining  the  position  of  com- 
puted points.  This  section  represents  the  modulus  of  the  function  along 
the  1 i ne  = 2.62. 


Qualitatively  Figure  6 shows  that  the  modulus  of  the  analytical 
quadratic  frequency  response  function  resembles  that  expected;  that  is, 
there  are  two  prominent  humps,  one  associated  with  difference  frequency 
i nteract i ons,  one  associated  with  sum  frequency  interactions,  and  in 
addition  the  two  humps  are  not  symmetrical.  At  bi-frequencies  distant 
from  those  of  the  main  humps  the  modulus  of  the  function  is  higher  than 


FIGURE  6 ISOMETRIC  VIEW  OF  THE  MODULUS  OF  THE  ANALYTICALLY  ESTIMATED 
QUADRATIC  FREQUENCY  RESPONSE  FUNCTION 
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was  expected  (though  it  must  be  admitted  that  the  expectation  was  not 
based  on  firm  experimental  estimates).  Clearly,  the  maximum  input  fre- 
quency should  have  been  extended  beyond  = 2.62  for  purposes  of  better 
defining  the  bi-frequency  domain  of  significance.  This  additional  com- 
putation was  not  carried  out  because  no  experimental  results  were 
available  for  the  corresponding  bi-f requenci es . 

The  experimental  estimates  as  presented  in  Ref.  1 are  not  quite 
equivalent  to  the  analytical  estimates  by  reason  of  a difference  in  phase 
reference  point.  The  phases  in  the  analytical  estimates  are  with  refer- 
ence to  wave  at  LCG.  The  phases  in  the  experimental  results  are  with 
respect  to  a wave  one  model  length  forward  of  LCG  which  had  been  time 
shifted  to  aid  in  the  cross-bi-spectral  analysis  process.  Within  the 
objectives  of  Ref.  1 it  was  convenient  to  retain  this  peculiar  phase 
reference.  For  present  purposes  it  was  convenient  to  correct  the  experi- 
mental results  to  the  same  phase  reference  as  had  been  assumed  for  the 
analytical  case.  The  net  phase  difference  between  the  time  shifted  wave 
and  wave  at  LCG  was  calculated  according  to  the  outline  in  Ref.  2.  This 
net  difference  is  equivalent  to  a linear  filter  with  unity  modulus  acting 
upon  the  input  to  the  system,  and  the  experimental  quadratic  frequency 
response  function  estimates  were  corrected  according  to  the  rules  in 
Refs.  1 and  5 for  cascaded  linear  and  quadratic  systems. 


The  qualitative  indications  mentioned  in  connection  with  Figure  6 
suggested  that  an  initial  comparison  between  the  magnitude  of  analytical 
and  experimental  estimates  might  best  becarried  out  indirectly.  The 
experimental  estimates  were  made  from  the  ratio  of  two  cross-bi-spectral 
estimates  in  accordance  with: 


A . \ A Pi 

Ga(<VU)a)  = Ga(— 2~ 


nA) 

o ' 


where  the  cap  denotes  an  estimate  and: 

) is  the  cross-bi-spectrum  of 
wave,  wave  and  resistance 

CTi(Q1,ns)  is  the  cross-b i -spectrum  of 
wave,  wave  and  (wave)2 . 


(52) 
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CTl^i^s)  “ ST|T1^U1i  ^ ST1T1^UJ2  ^ 

where 

S-q^(ai)  is  the  scalar  spectrum  of  wave, 

the  wave, wave,  (wave) 2 cross-rbi -spectrum  is  just  another  way  of  estimating 
the  product  of  scalar  wave  spectra.  When  C^(ni,Qa)  is  small  the  excita- 
tion of  the  system  at  that  bi-frequency  is  small,  and  the  estimates  implied 
by  Eq.  (52)  may  be  expected  to  be  seriously  corrupted  by  noise.  In  Refs. 

1 and  5 estimates  resulting  from  Eq.  (52)  where  C„(01,02)  was  less  than 
10%  of  its  peak  were  discarded.  Figure  7 is  an  isometric  plot  of  one  of 
the  experimental  wave, wave, (wave)2  cross-b i -spectra . This  cross-bi- 
spectrum (and  those  for  the  other  experimental  sea  states)  is  very  small 
in  much  of  the  bi-frequency  plane. 

These  considerations  suggested  that  a first  comparison  of  analysis 
and  experiment  should  be  made  by  computing  the  wave, wave, res i stance  cross- 
bi-spectrum  which  should  have  been  observed  by  experiment,  given  the 
analytical  results  for  the  quadratic  frequency  response  function  and  the 
experimental  wave  spectrum  or  the  experimental  results  for  C^(Q„ ,Q2) . The 
computation  for  the  modulus  of  the  wave, wave, res i stance  cross-bi-spectrum 
for  experimental  Sea  State  B for  instance,  amounts  to  multiplying  the 
function  shown  in  Figure  6 by  that  shown  in  Figure  7,  or: 

n,  +n  q_ 

| cDT(nx,Q2)  ! = [ g2T(-L2^  , -£2-l)  [ c^A,)  (53) 

in  which  the  subscripts  "T"  denote  analytical  estimates.  This  procedure 
has  the  effect  of  weighting  the  analytical  estimates  so  that  quadratic 
response  which  should  not  have  been  excited  in  the  experiments  does  not 
enter  in  the  comparisons. 

The  computation  defined  by  Eq.  (53)  was  carried  out  for  each  of  the 
three  experimental  sea-states.  These  "analytically"  estimated  cross-bi- 
spectra are  shown  along  with  the  corresponding  experimental  cross-b i-spectra 
in  Figures  8 through  10  in  isometric  views  having  the  same  conventions  as 
Figure  6.  At  the  top  of  each  figure  is  the  analytical  prediction. 
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The  triangular  bi-frequency  domain  corresponds  to  that  of  Figure  6. 

At  the  bottom  of  each  figure  is  the  modulus  of  the  experimental  cross- 
bi-spectrum.  Both  parts  of  each  figure  are  drawn  to  the  same  scale. 

On  the  whole,  in  Figures  8 through  10,  the  major  peaks  which  lie 

j 

near  the  axes  correlate  fairly  well.  The  largest  such  peaks  are  located 

at  approximately  the  same  sum  or  difference  frequencies  in  both  analyti- 

; cal  and  experimental  results.  The  difference  in  magnitude  is  about  30% 

for  Sea  States  A and  B and  for  the  peak  located  in  the  Q axis  in  Sea 

* 

j State  C.  The  peaks  along  the  sum  frequency  axis  for  Sea  State  C ( F i g— 

* ure  10)  are  different  by  about  a factor  of  2.  These  margins  of  difference 

are  roughly  the  same  as  that  shown  for  the  mean  resistance  operator  in 
F igure  4. 


A close  comparison  of  the  experimental  results  in  Figures  8 
through  10  discloses  a range  of  humps  lying  along  the  lineQ2  = 1.  The 
magnitude  decreases  as  the  sea  severity  increases.  This  characteristic 
is  not  present  in  the  analytically  derived  cross-b i -spectra . There  is 
in  addition  a significant  hump  along  the  lineOs  3.2  in  Figure  8 for 
Sea  State  A.  This  last  corresponds  to  the  contributions  of  a noise  fre- 
quency which,  as  indicated  in  Ref.  1,  has  little  to  do  with  added 
res i stance. 

Overall,  Figures  8 through  10  indicated  a fair  correlation  between 
analysis  and  experiment,  and  that  a more  detailed  comparison  of  real  and 
imaginary  parts  of  the  quadratic  frequency  response  functions  would  be 
worth-wh i 1 e. 

The  details  of  both  real  and  imaginary  parts  of  the  analytical 
quadratic  frequency  response  function  are  shown  in  comparison  with  the 
experimental  results  from  Refs.  1 and  5 in  Figures  11  through  16.  In  fig- 
ures 11  and  12  the  real  and  imaginary  parts  of  the  analytical  result  are 
compared  with  experimental  results  derived  from  model  Sea  State  A,  Figures 
13  and  14  contain  the  corresponding  comparison  with  results  from  model 
Sea  State  B,  and  Figures  15  and  16  involve  the  comparison  with  results 
from  Sea  State  C. 
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These  figures  are  essentially  plotted  tables.  The  magnitude  of 
the  function  involved  is  indicated  by  a number  plotted  in  the  bi-frequency 
plane.  The  analytical  results  are  shown  in  the  quadrant  bounded  by  the 
positive  Q1  and  Q„  axes.  In  order  to  show  the  experimental  results  in 
the  same  figure  the  octant  bounded  by  the  and  -tQ  axes  has  been  dis- 
placed in  the  negative  fi  direction.  The  experimental  results  shown  in 
this  octant  are  plotted  according  to  the  basic  symmetry  property  of  the 
quadratic  frequency  response.  If  analysis  and  experiment  were  in  exact 
agreement  the  experimental  results  would  be  an  exact  reflection  about  the 
Q2  axis  of  the  analytical  results  in  the  Q2  - octant.  A similar  pro- 
cedure was  adopted  to  show  the  experimental  difference  frequency  inter- 
actions. In  this  case  the  octant  bounded  by  the  and  -u^  axes  is 
displaced  in  the  negative  Qs  direction  and  the  experimental  values  are 
plotted  so  that  exact  ana  1 yt i ca 1 /exper imenta 1 agreement  would  require  the 
experimental  results  to  be  an  exact  reflection  about  the  axis  of  the 
analytical  results  in  the  octant.  (In  this  case  strict  compli- 

ance with  the  symmetry  properties  of  the  function  would  require  that 
the  sign  of  the  imaginary  part  of  the  experimental  result  be  changed. 

In  order  that  comparison  be  made  easier  this  has  not  been  done.) 

There  are  several  other  conventions  in  Figures  11  through  16. 

All  magnitudes  are  rounded  to  the  nearest  integer.  Analytical  results 
were  obtained  for  all  positions  in  the  - ni  quadrant  up  to  ^ = 2.62. 
However,  when  a numerical  result  was  between  -2.5  and  +2.5  it  was  omitted 

from  the  plot  — blank  spaces  in  the  analytical  result  may  be  taken  to 

mean  that  the  result  was  small.  There  are  two  possible  meanings  to 
blank  space  in  the  experimental  portion  of  the  figures.  The  first  is 
the  same  as  for  the  analytical  case.  The  second  is  responsible  for  the 
concentration  of  experimental  results  in  limited  portions  of  the  octants. 
It  is  that  all  experimental  estimates  which  resulted  from  a 
wave, wave, (wave)3  cross-b i -spectra  1 estimate  less  than  10%  of  its  peak 
value  were  discarded  as  being  of  dubious  statistical  merit.  (Actually 
there  were  indications  in  the  previous  work'*'’  that  this  truncation  level 
is  optimistic  for  the  sample  size  involved  — 25%  of  peak  could  well  have 

been  a better  choice.)  The  contours  shown  in  the  figures  are  for  the 

zero  1 evel . 
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The  results  may  be  conveniently  discussed  according  to  bands  of 
output  frequency,  £2  . The  first  band  of  frequency  to  be  considered  is 
that  above  QE  ^ 3.  In  this  band  the  magnitude  of  the  analytical  result 
generally  decreases  from  its  peak  along  the  axis  to  a minor  hump 
centered  at  approximately  = 4,  = 1.  Experimental  results  in  this 

band  are  shown  only  for  Sea  States  A and  B,  Figures  11  through  14,  and 
the  analyses  of  Ref.  1 indicated  that  much  of  the  content  of  this  band 
was  related  to  towing  system  induced  longitudinal  vibration  of  the  model. 
Thus  agreement  between  analysis  and  experiment  in  the  region  above  Qa  = 3 
was  not  expected.  Despite  this,  there  is  a surprising  resemblence  in 
the  position  of  the  zero  level  contours.  The  real  part  of  the  functions 
are  negative  at  Qs  = 3,  make  a transition  to  positive  halfway  between 
Q2  = 3 and  4 and  then  tend  to  return  to  negative  after  0 = 4.  The 

correspondence  in  zero  contour  position  in  the  imaginary  parts,  Figures 
12  and  14,  is  as  good,  a small  positive  area  appears  centered  on 
Qa  = 3.25  and  this  is  surrounded  by  negatives.  Agreement  in  magnitudes 
of  real  and  imaginary  parts  in  this  band  is  best  in  Figures  13  and  14 
for  Sea  State  B where  the  influence  of  vibration  on  experiments  was 
not  so  large.  The  corresponding  comparison  in  Figures  11  and  12  for 
Sea  State  A shows  large  differences  in  magnitude. 

The  fact  that  apparently  un-correlated  noise  appeared  to  influence 
the  cross-bi-spectral  estimates  of  response  could  not  be  explained  in 
Ref.  1.  The  similarity  in  the  behavior  of  the  analytical  and  experi- 
mental results  gives  rise  to  the  speculation  that  the  problem  vibration 
of  Refs.  1 and  5 might  have  been  induced  at  least  partially  by  the 
quadratic  component  of  the  wave  exciting  force  and  not  by  mechanical 
imperfections  in  rail  and  equipment  as  was  supposed. 

The  next  band  of  output  frequency  which  may  be  considered  is  that 
between  Qa  = 2 and  3.  This  is  the  band  where  the  experimental  estimates 
were  considered  most  significant.  The  level  of  agreement  between  the 
moduli  of  analytical  and  experimental  results  has  previously  been  noted 
in  connection  with  Figures  8-10.  The  position  of  zero  contours  in  the 
imaginary  part  of  all  experimental  estimates  correlates  quite  well  with 
the  analytical  (Figs. 12, 14  and  16) . However  in  the  real  parts  of  the 
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experimental  estimates  (Figs. 1 1, 13,  15)  the  zero  contour  is  consistently 
midway  between  the  limits  of  the  band  while  the  analytical  result  is 
almost  all  negative.  Here  there  is  a clear  inconsistency  between 
experiment  and  analysis.  The  fair  agreement  in  moduli  suggests  that 
the  phase  reference  correction  to  the  experimental  data  somehow  does  not 
reflect  the  actual  conditions  of  experiment,  or  that  the  radical  "hook" 
in  the  analytical  zero  contour  in  the  region  Qa  = 2,  = 1 should  not 

be  present.  Neither  supposition  can  be  proven  at  present. 

The  experimental  results  from  Sea  States  A and  B (Figs. 11-14) 
which  are  shown  for  Q2  less  than  2 in  the  upper  octant  correlate  badly 
with  the  analytical  result.  Only  in  the  case  of  the  imaginary  part 
derived  from  Sea  State  C (Fig. J 6)  is  there  reasonably  good  correlation. 

The  experimental  estimates  in  this  region  were  derived  from 
wave, wave, res i stance  cross-b i-spectra 1 estimates  of  relatively  low 
magnitude,  and  it  appears  that  the  analytical  result  bears  out  the  sup- 
position made  in  Ref.  1 that  the  experimental  results  in  this  region  tend 
to  suffer  from  sampling  variability. 

The  first  band  of  output  frequency  to  be  discussed  in  the  lower 
octant  is  that  from  fi  = 0 to  0.5  . For  practical  purposes  the  correla- 
tion of  real  parts  in  this  band  has  been  shown  in  Figure  4.  The  corre- 
lation of  imaginary  parts  in  this  band  may  be  regarded  as  fair,  considering 
that  the  imaginary  part  of  the  function  is  relatively  low  and  that  the 
experimental  results  are  sparse  in  some  cases. 

The  last  band  of  output  frequency  in  the  lower  octant  to  be  dis- 
cussed is  that  between  = 0.5  and  1.5.  This  is  an  interesting  band  in 
that  response  originating  in  quadratic  interactions  is  combined  with 
linear  response  at  frequencies  important  to  ship  motions.  The  analytical 
results  show  a significant  interaction  response  in  the  neighborhood  of 
= 1,  = 2.  The  experimental  results  from  Sea  State  B (Figs. 13, 14) 

come  closest  to  suggesting  the  same  interaction,  those  from  Sea  State  C 
(Figs. 15,16)  are  too  sparse  for  comparison,  but  those  from  Sea  State  A 
agree  rather  badly.  The  validity  of  some  of  these  experimental  results 
was  questioned  in  Ref.  1 since  they  also  were  derived  from  relatively 
low  values  of  wave, wave, res i stance  cross-b i -spectra. 
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In  the  remainder  of  the  = 0.5  to  1.5  band  a significant  response 
in  the  region  near  = 3 could  not  be  questioned  in  Ref.  1.  There  is 
clearly  no  correlation  with  the  analytical  result  in  this  case.  The  very 
significant  interaction  result  near  ^ = b which  was  obtained  by  reduct  ion 
of  data  from  Sea  State  A (Figs. 11, 12)  was  questioned  on  sampling  vari- 
ability grounds  in  Ref.  1,  and  it  is  clear  from  the  analytical  result  that 
this  interaction  cannot  be  substantiated. 
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SUMMARY 


The  general  objective  of  the  present  work  was  to  make  an  extension 
to  existing  hydrodynamic  theory  for  ship  resistance  added  by  waves,  to 
include  time  dependent  fluctuations  so  that  the  quadratic  frequency 
response  function  for  added  resistance  could  be  computed.  It  was  clear 
that  work  of  this  nature  might  provide  considerable  insight  into  the 
nature  of  resistance  fluctuations  in  random  seas,  and,  by  providing  a 
physical  model,  clarify  previously  obtained  experimental  estimates  of  the 
funct ion. 

Given  the  definition  and  interpretation  of  the  quadratic  frequency 
response  function  which  is  afforded  by  the  functional  polynomial  input- 
output  model,  as  well  as  an  existing  generalized  expression  for  hydro- 
dynamic  force  on  a body  in  potential  flow,  the  general  formulation  of 
the  hydrodynamic  problem  was  found  to  be  straightforward  within  the 
constraints  of  one  cardinal  assumption.  This  assumption  is  that  the 
dominant  contribution  to  second  order  non-linear  force  response  is  due  to 
the  products  of  first  order  potentials  which  arise  from  the  Bernoulli 
Equation  --  and  not  from  non  1 i near i t i es  in  the  potentials  themselves. 

The  resulting  general  formulation  is  not  limited  to  added  resistance, 
and  may  have  application  in  other  problems. 


The  evaluation  of  the  quadratic  frequency  response  function  for 
added  resistance  was  carried  out  according  to  a strip  method  in  which 
the  sectional  wave  diffraction  induced  forces  as  well  as  forces  induced 
by  body  motion  were  evaluated  by  close  fit  methods.  Because  the  mean 
added  resistance  operator,  which  is  generally  the  only  result  from  other 
analytical  methods,  is  a special  part  of  the  quadratic  response  function 
it  was  possible  to  compare  the  results  of  the  present  method  with  some 
other  analytical  results  and  with  results  from  the  usual  added  resistance 
experiment.  The  present  analytical  method  appears  valid  on  this  basis, 
and  in  fact  appears  to  have  promise  as  an  improvement  in  the  analytical 
computation  of  mean  added  resistance. 
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In  the  previous  experimental  work  three  sets  of  estimates  of  the 
quadratic  frequency  response  function  for  added  resistance  had  been  made 
for  each  of  two  model  speeds  of  a Series  60  0.60  block  model  by  means 

of  cross-b i -spectra  1 analysis  of  data  obtained  in  irregular  waves. 

Because  the  wave  severity  had  been  varied  over  a good  part  of  the  range 
of  practical  interest,  an  approximate  equivalence  of  the  main  features 
of  the  quadratic  frequency  response  functions  derived  from  data  obtained 
in  different  sea  states  was  taken  to  mean  that  the  concept  of  a quadratic 
frequency  response  function  was  viable.  There  were  necessarily  statis- 
tical complications  entering  into  the  data  reduction  process  so  that 
those  experimental  estimates  which  were  considered  to  be  at  least  tenta- 
tively valid  did  not  cover  the  entire  bi-frequency  domain  of  interest, 
and  were  not  themselves  all  of  the  same  reliability. 

The  present  analytical  methods  were  utilized  to  compute  the  com- 
plete function  for  one  of  the  conditions  of  the  experiments,  and  these 
results  were  compared  with  experiment  in  several  ways.  In  the  quantita- 
tive sense,  the  correlation  between  analytical  and  experimental  frequency 
response  functions  may  be  considered  fair,  at  least  for  the  predominant 
features.  In  this  context  it  should  be  remarked  that  the  correlation 
between  experimental  and  analytical  results  for  the  mean  added  resistance 
operator  can  also  only  be  considered  fair.  Better  quantitative  correla- 
tion for  the  rest  of  the  quadratic  frequency  response  function  can  hardly 
be  expected.  Over  all,  there  is  at  least  a reasonable  resemblence  between 
experiment  and  analysis  with  respect  to  the  contours  of  zero  level  in 
real  and  imaginary  parts  of  the  function;  that  is,  with  respect  to  phase. 
Analysis  and  experiment  do  not  agree  well  in  the  band  of  the  bi-frequency 
plane  in  which  the  interaction  response  is  at  frequencies  coinciding  with 
the  range  important  for  ship  motions.  Both  sets  of  results  imply  a sig- 
nificant interaction  response,  but  disagree  about  what  combination  of 
input  frequencies  are  of  most  importance. 

The  most  important  concern  in  the  present  work  was  qualitative. 

The  fundamental  questions  were:  a)  what  is  the  nature  and  general 

behavior  of  the  quadratic  frequency  response  function  outside  the  range 
of  experimental  data,  and  b)  could  the  predominant  features  of  the  func- 
tion as  shown  by  experiment  be  confirmed? 
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There  appears  to  be  good  qualitative  agreement  between  analysis 
and  experiment  with  respect  to  the  most  prominent  features  of  the  function. 
These  prominent  features  involve  a hump  containing  the  peak  values  of  the 
mean  added  resistance  operator,  and  a hump  of  typically  larger  magnitude 
associated  with  sum-frequency  and  second  harmonic  response.  Apart  from 
these  main  features,  the  analytical  results  imply  the  existence  of 
secondary  humps  in  many  of  those  regions  of  the  bi-frequency  plane  which 
could  not  be  covered  by  experiment.  Most  of  the  experimental  estimates 
which  were  previously  considered  of  marginal  reliability  were  also  not 
confirmed  by  the  analytical  results.  Since  it  was  these  marginal  experi- 
mental estimates  which  gave  rise  to  the  speculation  that  the  quadratic 
frequency  response  function  was  inordinately  complicated,  the  present  work 
indicates  that  the  function  is  merely  a moderately  complicated  function 
of  two  frequencies. 

Because  the  present  results  are  thought  to  be  the  first  of  their 
kind,  firm  conclusions  seem  premature.  It  is  apparent  that  confirmation 
of  the  actual  level  of  accuracy  of  the  analytical  method  must  await 
experimental  data  containing  fewer  caveats.  However  the  work  demonstrates 
that  it  is  feasible  to  make  an  at  least  approximately  correct  hydrodynamic 
estimate  of  the  response  function  --  which  in  turn  allows  interpretation 
or  prediction  of  non-linear  fluctuations  of  resistance  in  a random  seaway. 
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RECOMMENDATIONS 

It  appears  that  reasonable  qualitative  agreement  has  been  achieved 
between  experimental  and  analytical  estimates  of  quadratic  frequency 
response  for  added  resistance.  The  natural  next  step  might  be  to  refine 
the  computational  methods  so  as  to  achieve  better  detailed  agreement. 

In  the  present  case  however  the  only  experimental  estimates  which  were 
available  may  have  been  obtained  in  the  hardest  possible  way  --  that  is 
by  analysis  of  irregular  wave  experiments.  This  was  a natural  approach 
in  demonstrating  the  applicability  of  the  input-output  model,  but  the 
results  obtained  substantiate  the  earlier  fears  that  substantial  uncer- 
tanties  related  to  sample  size  exist.  It  is  accordingly  recommended 
that  the  most  fruitful  next  step  in  the  development  would  be  to  develop 
dual  frequency  experimental  techniques  (without  statistical  complications) 
which  would  produce  results  directly  comparable  to  analysis. 

The  analytical  methods  utilized  herein  appear  to  have  some  promise 
in  the  sense  of  improving  theoretical  estimates  of  mean  added  resistance 
operators.  Trial  application  of  the  method  to  a wider  range  of  ship 
forms  is  recommended. 
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APPENDIX  A 

DIFFRACTION  FORCES  AND  MOMENTS 
ON  A UNIFORMLY  TOWED  RESTRAINED  SHIP  IN  OBLIQUE  WAVES 

Q 

1 . A Brief  Description  of  the  Method  for  Zero  Speed 

The  problem  is  to  formulate  the  wave-exciting  forces  and  moments  on 
a fixed  ship  in  oblique  regular  waves. 

Set  the  right-handed  fixed  coordinate  system  X,Y,Z  on  the  ocean's 
surface.  Precisely,  the  X-Y  plane  rests  on  the  calm  water  surface  (see 
Figure  2).  An  incident  wave  which  progresses  across  the  X-axis  obliquely 
at  angle  p.  with  the  X-axis  is 


h = a e 


i (v  Xcosp  + v Ysinp-u)  t) 
o o o' 


where 


a = wave  ampl i tude 


i»>  = circular  wave  frequency 

v = wave  number  (=  u)2/g) 
o ' o 


The  velocity  potential  which  generates  the  wave  h is  expressed  in  the 


i ( v Xcosp-tu  t) 

*,(X,Y,ZS  Vo,fi;t)  * cp^Y.Zjv^e  0 


. ^ Z i v Y s i np. 

9,(Y,Z;vo^)  = - 121  e ° e (A-3) 

o 

The  diffraction  potential  is  regarded  as  the  direct  response  to  the 
incident  wave  potential.  In  view  of  this  fact,  we  assume  the  diffraction 
potential  to  be  similar  in  form  to  , as  follows: 


$D(X,Y,Z;volM,;t)  = cf0(Y,Z;vo>M.)e 


i (v  XCOSp-UD  t) 


where 
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1 


» 


I 


cp0(Y,Z;vo,n)  = J*  QOU)G(Y,Z;TU;vo)dC 
C(x) 


(A -5) 


In  the  above,  the  Green  function  G is  the  two-dimensional  pulsat- 
ing source  potential  of  unit  intensity  at  the  point  (T|,£)  in  the  lower 

9 

half  of  the  Y-Z  plane.  C(x)  indicates  that  the  integral  is  to  be  taken 
along  a hull-section  contour  below  the  calm  water  surface.  The  unknown 
source  intensities  Q.  which  are  discretely  distributed  along  the  con- 
tour C(x)  are  determined  by  satisfying  the  condition  of  non-penetration 
of  the  water  particles  across  the  contour.  This  kinematical  boundary  con- 
dition is 


dcpD  Scp( 

dn  dn 


on  C(x) 


(A -6) 


The  diffraction  potential  satisfies 

1.  the  continuity  condition 

2.  the  linearized  free  surface  condition 

3.  the  radiation  condition 

4.  the  deep  water  condition 


2 . Extension  of  the  Previous  Method  to  the  Non-Zero  Speed  Case 

In  dealing  with  the  diffraction  of  waves  due  to  a uniformly  moving 
ship,  it  is  required  to  introduce  the  moving  coordinate  system  x,y,z 
which  coincides  with  the  fixed  system  X,Y,Z  at  time  t=0  and  moves 
along  the  X-axis  at  speed  U (see  Figure  2).  The  transformation  of  the 
coordinate  system  is  then  given  by 

X = x + Ut 
Y = y 

Z = z (A-7) 

whereas  the  frequency  of  encounter  u)  is 

U)  = ^o  - U vocos^  (A-8) 

wi  th 

oj  = the  wave  frequency 

v = the  wave  number  uu2/q 
o o 


1 


Qi 


o: 
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The  incident  wave  and  its  corresponding  velocity  potential  in  the 
fixed  coordinate  system,  Eqs.(A-l.  A-2,  A-3)  , are  transformed  by  sub- 
stituting Eqs  .(A-7)  and  (A-8) • Thus, 


h = a e 


i (v^xcosp,  + v^ys  i np-uut) 


Mx,y,z;v  .njt)  =(f,(y,z;v>)e 


(v  XCOSp-UJt) 


v z i v ys i np 

cp1(y,z;vo,p)=-  e ° e ° (A- 1 1 ) 

o 

It  might  seem  that  the  diffraction  potential  which  is  regarded  as 
the  response  to  the  incident  wave  potent ial  ( A-10)  should  assume  an  ex- 
pression similar  to  that  given  in(A-«4).  However,  then  the  convective 

term  U f—  would  vanish.  The  time  derivative  of  the  diffraction 

ox  ot 

potential  with  respect  to  the  fixed  frame  is  now  replaced  by  the  sum  of 
local  and  convective  terms  ^ - U ^ with  respect  to  the  moving  frame. 
For  instance,  the  linearized  free  surface  condition  is 


vat  ■ u ZtJ  + 9 dTj  V*>y>z;  v.p.;t)  = o on  z=0 


In  view  of  the  above  remarks,  we  assume  that 


$0(x,y,z;  v ,p;t)  = cpD(x,y,z;  v ,p)e 


i ( v xcosp-uot) 


where 


90(x,y,z;v,p)  = cp“(y,z;v,p)  on  C(x) 


9d  = J*  Q.(Tl,C)G(y,z;Tl,C;v)dC  , 


v = ufVg 


* 
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and  cp“  varies  as  the  form  of  the  hull  section  varies  along  the  hull  length. 
The  above  assumption  paves  the  way  for  the  stripwise  calculation  to  be 
carried  out  in  a rational  manner. 

The  kinematical  boundary  condition  for  the  potential  cp^-Hp^  's 


dcpD(y>z;v.^)  d9|  (y.z;v0.M.) 


on  C(x) 


(A— 15) 


The  procedure  for  satisfying  the  kinematical  boundary  condition  follows. 
First,  define  the  odd  and  even  function  of  cp  in  Eq.(A-ll): 


/ , v z 

(o)  qa  o . , 

..  = ■2—  e s 1 n( 

uu  v 


voysinp,) 


(e)  ga  VoZ  . . . 

9,  = e cos(v  ysinp) 

o 

The  normal  derivatives  of  cpj°^  and  cp|6^  are  9'ven  by 


(A- 16) 


1 dtP|  ; V0Z 

— ^ = -id  e [ s i nps  i n( v^ys  i np)  s i ncrtcos ( v^ys  i np)  cosa] 


j dcp^  v z 

— uj^e  [ s i npcos  ( v^ys  i np)  s i no-s  i n( v^ys  i np)  coscr] 


(A-17) 


where  a is  the  slope  of  the  contour  C(x)  at  a point  (y,z) . 

On  substituting  (A-lA)  and  (A-17)  in  (A - 1 5 ) , we  have  for  the  sway- 
exciting force 

£ 0.-2M^+  £ Q^+j  = -a)Qe  ° l[sinpcos(vriy.sinp)sina. 


J 'J 


-sin(v  y.slnu)cosa. J 
v o 1 1 


-£  q(2>j(2>  + £ q<2>  if2)  =0 
j=i  J j=i  N+J  'J 


(A-l 8) 
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""  


- 


and  for  the  heave-exciting  force 


s ««>,<?>  ♦ s Q»>  j(?> 
j.i  j 'J  j-i  NiJ  ,J 


= 0 


- Zj  q(3)j(3)  + £ O^.3-  1 ^ = -u>  e 0 '[sinp,sin(v  y . s i np,)  s i na. 
j=1  J 'J  j=)  N+j  ij  O ^ V Oyi  ' 1 


+cos(v  y . s i nu,)  coso.  ] 
01  1 


(A-19) 


where  and  J^:?3  are  the  “influence  coefficients"  qiven  in  the 

ij  ij  s 


Appendix  of  Reference  9.  For  roll,  the  formula  is  the  same  as  for  sway, 
except  that  the  mode  number  is  4 instead  of  2. 


If  N sources  are  taken  (i.e.,  i=l,2,...  N) , we  obtain  a [2n]-[2nJ 
system  of  equations  with  2N  unknowns,  giving  the  real  and  imaginary  parts 
of  the  unknown  complex  source  intensities  Q.  . 


Solutions  of  ( A—  1 8)  and  (A-19)  yield  the  diffraction  potentials  9 


*(e) 


and  cpp  ' for  a strip  section  contour  C(x) 


The  d i f f ract i on- i nduced  hydrodynamic  pressure  at  a point  on  the  hull 
surface,  for  instance,  for  the  even  function  which  contributes  to  the 
heaving  force,  is 


p = jjiup  + pU  |^J  9^  (x,y,z;v,p,)e 


I V^xcosp 


(A-20) 


where  the  time  factor  e ,UJt  is  omitted. 


The  integration  of  the  vertical  components  of  the  above  pressure  first 
along  the  contour  C(x)  and  second  along  the  hull  length  between  -Xj  and  X, 
yields  the  d i ff ract ion- i nduced  heave-exciting  force  in  the  form 
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The  relation 


^7  I iP^6J  (x,y,z;v>M,)dy  = J*  i pu>  (x.y.z; v,pi)dy  (A-22) 

C(x)  C(x) 

has  been  used  in  the  foregoing  derivation. 

In  carrying  out  the  integral  along  the  contour  c(x)  in  Eq.(A-21),  we 
make  use  of  the  relation  (h-14)  for  cp"(y,z;v,n)  and  designate  the  integral 
by  fj(x)  , 


fS(x)  = J i puxpD^e^  (y)z;v,^)dy 
C(x) 

Now  the  integration  of  (A  — 21)  by  parts  with  respect  to  x yields 


D 

F,,  2 iv  xcosp.  _ „ -iv  A.cosu, 

J 6 ° ° 


(A- 23) 


(A -24) 


whe  re 


fS<-v 


is  Eq.(A-23)  for  x = 


The  di ffraction-induced  sway-  and  rol 1 -exc i t i ng  forces  and  moments 
are  evaluated  in  similar  fashion; 


D . 

Fl)h  r2  'V'”11  .0,  , . U ,D,  » , "Vl'”' 
-J--  J e fs<x)dx  "iW  fS('Ve 

“1 

D & 

f<ph  C2  ivoXCOS|1  n,  , . u ,D,  , , "'Vl'0*1* 

— ■ J,  e Vx)dx  - t:  V'V6 


where 


W = 


I*  . *( 

, ipuV 

:(x) 


fn(x)  = J ipuxpj^  (ydy  + zdz) 
C (x) 


(A -25) 


(A -26) 


(A-27) 


(A -28) 
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A similar  procedure  is  also  applied  to  calculate  the  pitch-  and  yaw- 
exciting moments.  For  example,  the  pitch-exciting  moment  is 


fd 

* h 
a 


COSLL  2 iv  XCOSLL  / \ 

xdx  ; iW<')d» 

-*] 


C(x) 


2 iv  xcosp.  . , , 

U r»  o , d p (e), 

- — e xdx  — ; 1 puxir  ' dy 

1 to  4.  dx  J . D 


-£ 


C(x) 


(A-29) 


On  substituting  f^x)  °f  Eq . (A  —23)  into  the  above  equation  and  integrating 
H 


by  parts  in  the  second  integral,  we  have 
.0 


Ffh 


X £ 

2 iv  xcosp,  n ,,  2 iv  xcosp,  . 

- j x d y*)d* + J e fH(x)dx 

-£]  -X, 


U£, 


uu  H'  r 


. -iv  X.cosp, 

D,  * n . o 1 


(A  — 30) 


The  yaw-exciting  moment  is  obtained  as 


.0 


Xh 

a 


x e 


IVqXCOS^ 


f^(x)  dx  - — f e 
S 1 (jo  J 

-l 


2 iv  xcosp.  _ 

0 *0,  \a 

f$(x) dx 


'1 


UX  -iv  £ cosu 

+ 1 f°(-Ve  ° ’ 


I UJ 


(A-31) 
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PRINCIPAL  NOTATION 


a.  amplitude  of  wave  component  of  frequency  uj. 


B 

beam  of  a section  in  the  load  waterplane 

C (x) 

a transverse  section  contour 

D (t) 
w 

resistance  added  by  waves 

— » 
F 

hydrodynamic  force 

F 

n 

Froude  Number 

fju,) 

sectional  heave  diffraction  force  per  unit 
wave  ampl i tude 

U) 

linear  frequency  response  function 

quadratic  frequency  response  function 

gx(t) 

linear  impulse  response 

second  degree  kernel  analogous  to  an  impulse  response 

\ (<JJj  ,U3j  ) , 

H!  U;>-Wj) 

functions  defined  by  Eqs . (20)  and  (21) 

H2 

H2  ( a;  j j -auj  ) 

functions  defined  by  Eqs.  (22)  and  (23) 

L 

hull  1 ength 

mM(u>.) 

sectional  heave  added  mass 

N(U3.) 

sectional  damping  coefficient 

n 

outward  unit  normal  vector 

n 

the  z-component  of  the  unit  normal  on  the  hull  surface 

— * 
R 

quadratic  part  of  the  added  resistance 
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force  due  to  two  superimposed  regular  wave  trains 

submerged  surface  of  the  body 

portion  of  the  free  surface  inside  a far  field 
control  surface,  S 

CO 

average  section  draft 


ship  speed 


v(m)(«.) 


velocity  amplitude  of  heave  (m=3)  and  pitch  (m=5) 


X,Y,Z 


fixed  coordinates 


x,y,z 


moving  coordinates 


model  displacement 


heave  ampl i tude 


wave  elevation  of  the  i wave  train 


wave  elevation 


wave  number  corresponding  to  i wave  train 


mass  density 


total  velocity  potential 


body  potential 


incident  wave  potential 


wave  resistance  potential,  steady  state 


time  dependent  potential  due  to  waves 
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CPB(<\)  Exp[-iu)kt] 

CPD(u,i) 

cpj  (u)k)  Exp[-  'Vl 

♦ 

q2 

,<u2 


body  potential,  encounter  frequency  uj^ 

the  diffraction  potential  per  unit  wave  amplitude 
in  the  i wave  train 

incident  wave  potential,  encounter  frequency,  u>k 

radiation  potential  due  to  heaving  and  pitching 
motions,  each  per  unit  velocity  amplitude 

pi tch  ampl i tude 

difference  frequency 

sum  frequency 

frequency 

t h 

wave  frequency  corresponding  to  i c wave  train 
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1 McDonnell  Douglas  Aircraft  Co. 

3855  Lakewood  Blvd. 

Long  Beach,  CA  90801 
Attn:  T.  Cebeci 


1 U.S.  Merchant  Marine  Academy 
Kings  Point,  L.  I .,  NY  11204 
Attn:  Academy  Library 

1 Bolt,  Beranek  „ Newman 

50  Mou 1 ton  S+  reet 
Cambridge,  MA  02138 
Attn:  Library 

1 Bethlehem  Steel  Corporation 

Center  Technical  Division 
Sparrows  Point  Yard 
Sparrows  Point,  MD  21219 

1 Bethlehem  Steel  Corporation 

25  Broadway 
New  York,  NY  10004 
Attn:  Library  (Shipbuilding) 

» 


1 Lockheed  Missiles  & Space  Co. 
P.0.  Box  504 
Sunnyvale,  CA  94088 
Attn:  Mr.  R.L.  Waid,  Dept  57-74 

Bldg.  1 50,  Fac i 1 i ty  1 


1 Newport  News  Shipbuilding  & 
Dry  Dock  Company 
4101  Washington  Avenue 
Newport  News,  VA  23607 
Attn:  Technical  Library  Dept. 


1 Nielsen  Engineering  6-  Research  Inc 

510  Clude  Avenue 
Mountain  View,  CA  94043 
Attn:  Mr.  S.  Spangler 

1 Ocean ics,  Inc. 

Technical  Industrial  Park 
P la  inview,  L.I.,  NY  11803 


» 
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Society  of  Naval  Architects 
and  Marine  Engineers 
74  Trinity  Place 
New  York,  NY  10006 
Attn:  Technical  Library 

Sun  Shipbuilding  & Dry  Dock  Co. 

Chester,  PA  19000 

Attn:  Chief  Naval  Architect 

Sperry  Systems  Management 
Division 

Sperry  Rand  Corporation 
Great  Neck,  NY  11020 
Attn:  Technical  Library 

Stanford  Research  Institute 
Menlo  Park,  CA  94025 
Attn:  Library  G-021 

Southwest  Research  Institute 
P.0.  Drawer  28510 
San  Antonio,  TX  782 8k 

Attn:  Applied  Mechanics  Review 

Dr.  H.  Abramson 

Tracor,  Inc. 

6500  T racor  Lane 
Austin,  TX  78721 

Mr.  Robert  Taggart 
3930  Walnut  Street 
Fairfax,  V A 22030 

Ocean  Engr.  Department 
Woods  Hole  Oceanographic  Inst. 
Woods  Hole,  MA  02543 

Worcester  Polytechnic  Inst. 
Alden  Research  Laboratories 
Worcester,  MA  01609 

Applied  Physics  Laboratory 
University  of  Washington 
1013  N.E.  40th  Street 
Seattle,  WA  98105 
Attn:  Technical  Library 


University  of  Bridgeport 
Bridgeport,  CT  06602 
Attn:  Dr.  E.  Uram 

University  of  California 
Naval  Architecture  Department 
College  of  Engineering 
Berkeley,  CA  94720 
Attn:  Library 

Prof.  W.  Webster 
Prof . J . Pau 1 1 i ng 
Prof.  J.  Wehausen 

California  Institute  of  Technology 
Pasadena,  CA  91 1 09 
Attn:  Aeronautics  Library 

Dr.  T.Y . Wu 
Dr.  A.J.  Acosta 

Docs/Repts/Trans  Section 
Scripps  Institution  of 
Oceanography  Library 
University  of  California, 

San  Diego 

P.0.  Box  2367 

La  Jolla,  CA  92037 

Catholic  University  of  America 
Washington,  DC  20017 
Attn:  Dr.  S.  Heller,  Dept  of 

Civil  & Mech  Engr 

Florida  Atlantic  University 
Ocean  Engineering  Department 
Boca  Raton,  FL  33432 
Attn:  Technical  Library 

Dr.  S.  Dunne 


Univers i ty  of  Hawai i 
Dept,  of  Ocean  Engineering 
2565  The  Mall 
Honolulu,  HI  96822 
Attn:  Dr.  C.  Bretschneider 


Institute  of  Hydraulic  Research 
The  University  of  Iowa 
Iowa  City,  I A 52240 
Attn:  Library 
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Department  of  Ocean  Engineering  2 

Massachusetts  Insti.  of  Technology 
Cambridge,  MA  02139 
Attn:  Department  Library 

Prof.  P.  Mancie! 

Prof.  M.  Abkowitz  , 

Dr.  J.  Newman 

St.  Anthony  Falls  Hydraulic 
Laboratory 

University  of  Minnesota 
Mississippi  River  at  3rd  Ave.  S.E. 
Minneapolis,  MN  55414 
Attn:  Prof.  E.  Silberman  ^ 

Dr.  C.  Song 

Department  of  Naval  Architecture 
and  Marine  Engineering 
University  of  Michigan 
Ann  Arbor,  Ml  48104 
Attn:  Library  ^ 

Dr.  T .F . Og i 1 v i e 
Prof.  F.  Hammitt 

College  of  Engineering 

University  of  Notre  Dame  . 

Notre  Dame,  IN  46555 

Attn:  Engineering  Library 


Stanford  University 
Stanford,  CA  94305 
Attn:  Engineering  Library 

Dr.  R.  Street 

Webb  Institute  of  Naval 
Archi tecture 
Crescent  Beach  Road 
Glen  Cove,  L. 1 .,  NY  11542 
Attn:  Library 

Prof.  E.V.  Lewis 
Prof.  L.W.  Ward 


Applied  Research  Laboratory 
P.0.  Box  30 

State  College,  PA  16801 
Attn:  Dr.  B.  Parkin,  Director 

Garfield  Thomas 
Water  Tunnel 

Dr.  Michael  E.  McCormick 
Naval  Systems  Engineering  Dept. 

U.S.  Naval  Academy 
Annapol is,  MD  21402 

Dr.  Douglas  E.  Humphreys  (Code  712) 
Naval  Coastal  Systems  Laboratory 
Panama  City,  FL  32401 


Davidson  Laboratory 

Stevens  Institute  of  Technology 

71 1 Hudson  Street 

Hoboken,  NJ  07030 

Attn:  Library 

Dr . J . Bres 1 i n 
Dr.  S.  Tsakonas 


